Anisotropies in thermal Casimir interactions: ellipsoidal colloids 

trapped at a fluid interface 

Ehsan Noruzifar and Martin Oettel 
Institut fur Physik, WA 331, Johannes- Gutenberg- Universitat Mainz, 55099 Mainz, Germany 

(Dated: April 21, 2009) 

Abstract 

We study the effective interaction between two ellipsoidal particles at the interface of two fluid phases 
which are mediated by thermal fluctuations of the interface. In this system the restriction of the long- 
ranged interface fluctuations by particles gives rise to fluctuation-induced forces which are equivalent 
to interactions of Casimir type and which are anisotropic in the interface plane. Since the position 
and the orientation of the colloids with respect to the interface normal may also fluctuate, this system 
is an example for the Casimir effect with fluctuating boundary conditions. In the approach taken 
here, the Casimir interaction is rewritten as the interaction between fluctuating multipole moments 
of an auxiliary charge density-like field defined on the area enclosed by the contact lines. These 
fluctuations are coupled to fluctuations of multipole moments of the contact line position (due to the 
possible position and orientational fluctuations of the colloids). We obtain explicit expressions for the 
behavior of the Casimir interaction at large distances for arbitrary ellipsoid aspect ratios. If colloid 
fluctuations are suppressed, the Casimir interaction at large distances is isotropic, attractive and long 
ranged (double-logarithmic in the distance). If, however, colloid fluctuations are included, the Casimir 
interaction at large distances changes to a power law in the inverse distance and becomes anisotropic. 
The leading power is 4 if only vertical fluctuations of the colloid center are allowed, and it becomes 8 
if also orientational fluctuations are included. 

PACS numbers: 82.70.Dd, 68.03.Kn 
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I. INTRODUCTION 



When a fluctuating medium with long-ranged, power-law correlations is confined between a 
set of boundaries, forces with likewise long-ranged character are induced between the bound- 
aries. There are different possible sources of such fluctuations in a medium: in a quantum- 
mechanical system it is the zero point energy of the vacuum (or ground state), and in a classical 
system it is the finite temperature which causes order parameter fluctuations [l]. 

This kind of force was discovered theoretically by Casimir in 1948 for the case of two par- 
allel, conductive and uncharged plates immersed in vacuum which he attributed to zero-point 
fluctuations of the electromagnetic field For a review on recent progress and the status of 
experimental verification of this quantum mechanical Casimir effect, see Ref. A classical 
equivalent of the Casimir force observable between objects immersed in a fluid in the vicinity of 
its critical point was predicted 30 years later |4|. The fluctuations of the order parameter field 
near the critical point are long-ranged, and thus they give rise to a Casimir-like, fluctuation- 
induced force. This effect has recently been observed in an experiment probing the force on 
colloidal particles which have been immersed in a near-critical binary mixture in the vicinity of 
a wall Q|. Another classical variant of the Casimir interaction is found between particles (col- 
loids) that are trapped at membranes 



. In this 



7| or at the interface of two fluid phases 
two-dimensional latter instance, thermally excited height fluctuations of the interface which 
have long-ranged nature are disturbed by the presence of colloids. The energy spectrum of the 
fluid interface on a coarse-grained level is very well described by an effective capillary wave 
Hamiltonian which governs both the equilibrium interface configuration and the thermal fluc- 
tuations around this equilibrium. Since capillary waves are the Goldstone modes of the broken 
translational symmetry pertaining to a free interface, their correlations decay logarithmically 
in the absence of gravity and the corresponding fluctuation-induced forces are a manifestation 
of the Casimir effect for a Gaussian theory in two dimensions. Compared to manifestations of 
the Casimir effect in a three-dimensional bulk medium, a new phenomenon arises here: The 
boundary of the fluctuating interface, which is represented by the contact line on the colloid 
surface, is itself mobile due to position fluctuations of the colloids and finite surface tensions 
of the colloid-liquid interfaces and thus the Casimir force receives another contribution due to 
these fluctuating boundaries. This effect has been noticed first for colloidal rods trapped at 
membranes and films For a system composed of two spherical colloids trapped at a fluid 
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interface, various realizations of these fluctuating boundary conditions are possible and it has 



been shown 



111 ] that the fluctuation-induced force sensitively depends on the type of 
boundary conditions, with the asymptotics of the force ranging from l/(d\nd) to d~ 9 . 

In recent work, the general form of the fluctuation- induced interactions between a finite num- 
ber of compact objects of arbitrary shapes and separations has been calculated for a fluctuating 
medium of scalar Gaussian type 12| and an electromagnetic medium 13| with fixed boundary 
conditions on the surface of the objects. This has been achieved by viewing the Casimir interac- 
tion as resulting from fluctuations of source distributions (of the fluctuating field) on the surfaces 
of the objects which are decomposed in terms of multipoles. Then by a functional integral over 
the effective action of multipoles, the resulting interaction has been found. In such a way, the 
effect of anisotropic object shape on the Casimir force appears to be tractable. However, the 
objects' shape and the boundary conditions enter the effective action by its scattering matrix 
which is a nontrivial object already for simple shapes 12|]. Studies with explicit calculations 
for objects other than spheres, cylinders and walls have partially been focused on the effect 

n n n □ 

of wall corrugations |14l . 1151 ] , but also sharp edges pi| and rectangular bodies [17J have been 
investigated. Also, in a recent work, the quantum Casimir interaction between two ellipsoidal 
particles as well as an enipsoida, p artlC ,e and a plane has been stndied Q. In the present 
work, we investigate the fluctuation-induced interactions between two ellipsoidal colloids that 
are trapped at the interface of two fluid phases with special emphasis on the effect of anisotropy. 
Ellipsoids (spheroids) of varying aspect ration allow a smooth interpolation between spheres [9( 
and rods [3]. The ellipsoidal colloids are assumed to be of Janus type, therefore the interface 
contact line is always pinned to the colloids surface, nevertheless the vertical position of colloids 
and their orientation may fluctuate, giving rise to the already mentioned feature of fluctuating 
boundaries. For the calculation of the Casimir force, we employ techniques which have been 
introduced in previous work and which partially can also be interpreted in terms of 



the scattering matrix ansatz of Refs 



13| such that our results concrete example of 



the general theory for the effects of anisotropic object shape. As stressed before, however, at 
an interface we always have the influence of the fluctuating boundary conditions. In order to 
study the effect of the mobile boundaries in our work, we have divided our investigations into 
two main parts; an interface fluctuation part and a second part where the colloid fluctuations 
are included. In the interface fluctuation part, the positions of the colloids and thus the contact 
lines are fixed and therefore the problem reduces to the "usual" Casimir problem with Dirichlet 
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boundary conditions. The Casimir interaction is determined by multipoles of an auxiliary field 
on the (Id) contact lines which is similar to the source field in the language of Ref. [l^ . In the 
second part, we include that the colloid position may fluctuate in all possible ways (height and 
tilts) which turns out to lead to Casimir interactions determined by multipoles of an auxiliary 
field defined on the 2d domain enclosed by the contact lines. The fluctuating boundaries result 
in certain restrictions on these multipoles. (Incidentally, we note that for small fluctuations 
amplitudes, an exact separation of interface and colloid fluctuation in the partition sum is pos- 
sible (9), [lo|. This route, however, appears to be very difficult to follow for other than spherical 
colloids and is not taken here.) 

The paper is structured as follows. In Sec. [TT1 we introduce the effective Hamiltonian of the 
model system and the partition function which is defined by functional integrals over colloid and 
interface fluctuations. The implementation of the different boundary conditions is discussed. In 
Sees. Ill II and IIII A~| we determine the fluctuation-induced force for the interface fluctuation part 
in the long and short distance regimes analytically and for intermediate distances numerically, 
respectively . In Sec. IIVI we include the colloid fluctuations to the problem and discuss the 
modified long-distance asymptotics of the Casimir interaction. Sec.|V]contains a brief summary. 
Some technical details of the calculations have been relegated to Apps. lATlDl 



II. MODEL 



The investigated system consists of two nano- or microscopic, uncharged spheroidal colloids 
with principal axes a, b, b (a > b), which are trapped at the interface of two fluid phases I and 
II. The effective interaction between the colloids is mediated by thermal height fluctuations of 
the (sharp) interface. Without fluctuations, the equilibrium interface is flat and is set to be at 
z = 0. The corresponding equilibrium position of the colloids is assumed to be symmetrical with 
respect to z — > —z, such that at the contact line the contact angle is ir/2. The elliptic cross- 
section of the ellipsoids with the equilibrium interface is denoted by S^rcf which is an ellipse 
with major and minor axes a, b, respectively. SVref may also be expressed in confocal elliptic 
coordinates by the elliptic radius £o> see App. |X]for the coordinate definitions. The equilibrium 
interface at z = without the two elliptic holes SVref cut out by the colloids is termed the 
reference meniscus S'mcn.ref = R 2 \ UjfS^ref. Deviations from this planar reference meniscus are 
considered to be small, without overhangs and bubbles, therefore the Monge representation 
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FIG. 1: Top view on the system 



(x,y,z = u(x,y)) = (x, z — u(x)) is employed to describe the interface position. The colloids 
are of Janus type, thus the contact line is always pinned to their surface. The total Hamiltonian 
of the system which is used for calculating the free energy costs of thermal fluctuations around 
the flat interface is determined by the change in interfacial energy of the interface I /II: 



Where 7Av4 mcn expresses the energy needed for creating the additional meniscus area associated 
with the height fluctuations. In Eq. (CQ), S mcn is the meniscus area projected onto the plane 
z = (where the reference interface is located) and .Smen.ref is the meniscus in the reference 
configuration mentioned above. The first line of Eq. ([!]) constitutes the drumhead model which 
is well-known in the renormalization group analysis of interface problems, but is also used for 
the description of elastic surfaces (c.f. Ref. In the second line we have applied a small 

gradient expansion which is valid for slopes |Vit| <C 1 and which provides the long wavelength 
description of the interface fluctuations we are interested in. The small gradient expansion 
entails that 



is the change in projected meniscus area with respect to the reference configuration. We rewrite 
this change in projected meniscus area in terms of the interface position f\ = u(dSi^ TC {) at the 
reference contact line ellipses dSi tVe f. fi corresponds (in second order approximation) to the 
contact line of the colloid i with fluctuating center position hi and fluctuating orientation. The 
contact line which is a function of the elliptic angle t] only (see App. |A] for the definition of 
elliptic coordinates) is expanded as 




(1) 






(3) 



m=0 
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and we refer to the coefficients Pj m and Qi m as boundary multipole moments below. The 
desired expression of AAproj in terms of boundary multipole moments proceeds as discussed in 
Ref. and allows us to identify it as a sum over boundary Hamiltonians Ti.^ for each colloid 
i (see also App. IB]) : 

7 AA proj = £W b ,i[/i] 

i 

= J2y ( tanh ^ p l + coth ^ <&) • ( 4 ) 

i 

Putting Eqs. (CQ) and (j4j) together, the total change in interfacial energy is the sum 

H tot = H cw + H htl + H h , 2 = | / d 2 x (Vw) 2 + W b ,i + ^b, 2 (5) 

^men.ref 



of the capillary wave Hamiltonian 7^ cw which describes the energy differences associated with 
the additional interfacial area over the reference configuration and the boundary Hamiltonians 
7Yb,i which can be viewed as the energy cost due to fluctuations of the contact line (and which 
in turn are caused by colloid height and tilt fluctuations). As is well-known, the Hamiltonian 
7i cw is plagued with both a short- wavelength and a long-wavelength divergence which, however, 
can be treated by physical cutoffs. The natural short- wavelength cut-off is set by the molecular 
length-scale a of the fluid at which the capillary wave model ceases to remain valid. The long 
wavelength divergence is reminiscent to the fact that the capillary waves are Goldstone modes. 
Of course, in real systems the gravitational field provides a natural damping for capillary 
waves. Accounting also for the costs in gravitational energy associated with the interface 
height fluctuations, therefore, introduces a long wavelength cutoff and leads to an additional 
term ("mass term") in the capillary wave Hamiltonian, 

,2" 

(6) 



2 



= Z I d 2 x 



(v«) 2 + £ 



'-'men, ref 

Here the capillary length is given by A c = [7/dpn — Pilg)} 1 ^ 2 , where is the mass density 
in phase i and g is the gravitational constant. Usually, in simple liquids, A c is in the range 
of millimeters and, therefore, is by far the longest length scale in the system. In fact, here 
it plays the role of a long wavelength cutoff of the capillary wave Hamiltonian H. cw , and we 
will discuss our results in the limit A c ^> R and A c ^> d. However, as we will see below, 
care is required when taking the limit A c — ► 00 (corresponding to g — > 0), since logarithmic 



divergencies appear 



211 ] . Another common way to introduce a long- wavelength cut-off is the 
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FIG. 2: Side view on the system 

finite size L of any real system. As discussed in Ref. [ijj], the precise way of incorporating 
the long-wavelength cut-off is unimportant for the effects on the colloidal length scale. As 
an example, in both approaches the width of the interface related to the capillary wave is 
logarithmically divergent, (w(0) 2 ) ~ lnA c [L]/cr. Via the integration domain of T~C CW , the total 
Hamiltonian of the system, Eq. ([5]), implicitly depends on the geometric configuration. This 
leads to a free energy T{d, 6\, 6 2 ) which depends on the distance d between the colloid centers 
and the orientation angles Q\ and 9 2 of their major axes with respect to the distance vector 
joining the colloid centers (see Fig. [1]). The free energy is related to the partition function 
Z(d, 61,62) of the system by 

F(d, #1, 6 2 ) = ~k B T In Z(d, 0i, 62) (7) 

The partition function is obtained by a functional integral over all possible interface configura- 
tions u and boundary configurations ff, the relation between interface and boundary configu- 
rations is included by ^-function constraints, 

Z = Z^ J Vu exp {-^i} fl J V f* II 5 Hx.) - /,(*)] exp 

i=l X 4 e9Si,ref 

Here Zq is a normalization factor such that Z(d 00) = 1 and ensures a proper regularization 
of the functional integral. Via the (^-functions the interface field u is coupled to the contact 
line height fa and therefore, the boundary Hamiltonians Ti.^ have a crucial influence on the 
resulting effective interaction between the colloids. 

The kind of possible contact line fluctuations ft is solely determined by the colloid fluctuations 
since the contact line is pinned. These fluctuations are vertical fluctuation of the colloids on the 
axis normal to the equilibrium interface (height) and orientational fluctuations around that axis 
(tilts). In order to incorporate various boundary counditions into the solutions, we categorize 
them into three cases: 

(A) colloids are fixed in the reference configuration, thus there are no integrations over the 
boundary terms. 
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(B) colloid heights fluctuate freely without tilting, thus the boundary monopoles must be 
included in the integration measure so that T>fi = dPiQ. 

(C) unconstrained height and tilt fluctuations. Up to second order in the tilts this corresponds 
to the inclusion of boundary dipoles in the integration measure, thus T>fi = dP^ dPn dQn. 

Case (A) corresponds to the "standard" Casimir effect in 2d with Dirichlet boundary conditions 



fi = u(dSi tTC f) = 0. We call this the interface fluctuation part and it will 
Note that a short summary of this part has already been given in Ref. 
the colloid height and tilt fluctuations in (B) and (C) is given in Sec. IIVI 

III. INTERFACE FLUCTUATION PART 



De treated in Sec. IIHI 



241 ]. The inclusion of 



The partition function Z m for fixed contact lines fi — is given by 

(9) 



Z- m = Zv 1 Jvuf[ Y[ 6(u(xi)) exp | - 

i=l TUedSirei 



He™ \u] 



k B T 



The disapearance of the interface fluctuations at the colloids boundaries is included by the 
Dirac delta function. In this section, analytical expressions for the fluctuation induced force in 
the intermediate asymptotic regime a <C d <C A c are calculated. We express the 5-functions 
in Eq. by their integral representation via auxiliary fields ^j(xj) defined on the reference 
contact lines dS ijTe {. This enables us to integrate out the field u leading to 

Z, n = Z'- 1 [ rfp^expl-^ V / d£i [ di^^Gd^-^l)^)] ,(10) 

J i=l { 27 iJ=l JdS i,ret JdSj.r* J 

where dli is the infinitesimal line segment on the circles dSi tie f. After this integration, the 
normalization factor is changed, Zq — > Zq, such that still Z- m (d — > oo) = 1 holds. In Eq. (fTUl) 
we introduced the Greens function of the operator (—A + A^ 2 ) which is given by G(x) = 
Ao(|x|/A c )/(27r) where Kq is the modified Bessel function of the second kind. In the range 
d/\ c «C 1 and r /A c 1, we can use the asymptotic form of the K for small arguments, 
such that 27rG(|x|) £s — ln(7 e |x|/ 2A C ). Here, 7 e w 1.781972 is the Euler-Mascheroni constant 
exponentiated. We introduce auxiliary multipole moments as the Fourier-transforms of the 
auxiliary fields ipi on the reference contact line dSi, ie {, 

^ r 2n 

$im = dr) i h(r] i )cos(mT] i )il; i (-}c i (T] i )) , 



itm = / *7i sin(m^) Vi(xi(^)) , (11) 
Jo 

where r/j is the elliptic angle pertaining to a coordinate system centered around each colloid 
i, respectively, such that the x-axis in this colloid-specific coordinate system joins the two 
foci of S^ref. Furthermore, h(r)i) is the scale factor in elliptic coordinates (see App. |A]). The 
lengthy calculation leading to the multipole (Fourier) decomposition for the Greens function 
G(|xj — Xj|) (for general orientations 9\ and 82 of the ellipsoids) is given in App. O The final 
results is collected in Eq. fICll) and Eqs. flC15j) - f)C17j) . Using this, the double integral in the 



exponent of Eq. (fTOj) can be written as a double sum over the auxiliary multipole moments 
(Fourier components), consisting of a self-energy part G S cif when Xi and Xj reside on one ellipse 
and G- m t when the points X{ and Xj reside on different ellipses, respectively. The functional 
integral over the auxiliary fields becomes a product of integrals over their multipole moments, 
Vipi = dipio YlJLi dipijdijjfj, and the resulting partition function then reads 




where the vectors SE^ = V&f) with \I>^ = (-?/^ , 7/^, . . .) and = ■ ■ •) contain the 

auxiliary multipole moments of colloid i. The coupling matrix G which contains the Fourier 
modes of the Greens function G(xi — Xj) has a block structure. The self energy submatrix G se if 
which describes the coupling between auxiliary moments of the same colloid are diagonal, and 
its form can be determined from definition (fTTl) and Eq. fICll) . 

27r(^) T G sclf ^ = -ln^— (^ ) 2 + 2^ cosh«o)(^J 2 + sinh« )(^J 2 , (13) 

c n=l 

where a' 2 = a 2 — b 2 . The offdiagonal blocks Gj nt characterise the interaction between the 
multipole moments residing on different colloids. It is convenient to split the matrix into a 
block structure describing the interaction of cosine and sine multipoles: 

2vr (* X ) T G int § 2 = I I I I Z I (14) 



ss 
' int 




The matrix elements of the such defined sub matrices follow from Eqs. (IC15I) - (IC17I) . and are 
explicitly given by: 

G£L = -^(^f] (15) 



00 V 2A 



mn 

1=0 



i \ rn+n+2l 

a 



E ( Td I A<mnl{61 ' ° 2) cosh ^ cosh «o) (16) 



(l \ m+n+2l 
^ A c mnl (9i, 9 2 ) smh(m&) sinh« ) (17) 



(I \ m+n+zl 
^ A^JB U 9 2 ) sinh(meo) cosh« ) (18) 

From Eq. (|T2j) we find that the fluctuation part of the free energy reads 

Fm = —k B T In Z in = — ln(det G) + const. , (19) 

,c[s] 



where const. = — kBTlnZo'. The factors #2), given in Eqs. (IC16j) and (IC17p . contain 



the dependence on the orientation angles $i and #2 of the ellipsoids (see Fig. [T]). As can be 
seen from above, the interaction coefficients (G^ c j between multipoles of order m and 

V / mn 

n take the form of a series in starting at l/d m+n . (For spherical colloids, this multipole 
interaction coefficient only contains the order l/d m+n [lo|.) In principle, the matrix G is infinite 
dimensional and det(G) is divergent and its regularisation is provided by the normalization 
factor Z' Q . The explicit series for the elements of G; nt allows for a systematic expansion of the 
logarithm in Eq. (fl9l) in powers of 1/d, 

2 

A 

where the coefficients j l 2 \ depend on the logarithms — ln(7 e d/ 2A C )) and — ln(7 e a'e^ )/ 8A C ), as 
well as the angles Q\ and 9 2 . The number of auxiliary multipoles included in the calculation of 
the asymptotic form of ^-fl uc in Eq. ( J20l) is determined by the desired order in 1/d. Inclusion of 
multipoles up to order n leads to an asymptotics correct up to l/d 2n . In the limit X c /d — > 00 
the free energy expansion coefficients in Eq. (120]) up to fourth order are 1 



= k E Tj2fZ (V) , (20) 



fin 

J2 ~ 



2 In 



Ad 



a 



{a + bf , 3 , 



16 32 



(a 2 -6 2 )(cos(20i) +cos(20 2 )) 



(21) 



1 Note that a factor of l/(21n[2<i/(a + b)]) is missing in the expression for f^ 1 in Eq. (29) of Ref. [24 1. 
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fill 

2 11 



1 



In 



Ad 



16(a - 6) (a + 6) 3 (cos(20i) + cos(2# 2 )) 



,a + 6, 

+ ll(a 2 - 6 2 ) 2 (cos(4^) + cos(4# 2 )) 
+44(a 2 - 6 2 ) 2 cos(20i + 2fl 2 ) + 6(a + 6) 4 

+ [8(a 2 - 6 2 ) 2 cos(2^ + 20 2 ) + 8(a + 6) 4 ] 



The double-logarithmic divergence in d in the leading coefficent f$ n is a reflection of the fact 
that the interface itself becomes ill-defined for A c — > oo due to the capillary waves. For the 
Casimir force itself, however, we find a finite value for all d in the limit A c — > oo. Anisotropies in 
the Casimir interaction appear here first in the subleading term / 2 n . Their angular dependence 
stems from the monopole-dipole interaction of the auxiliary field, and the attraction is maximal 
if both ellipses are aligned tip-to-tip. 

In the opposite limit of small surface-to-surface distance h — d — da <C dd, (where d c \ is 
the distance of the closest approach between ellipses) the fluctuation force can be calculated by 
using the Derjaguin (or proximity) approximation 22j. It consists in replacing the local force 
density on the contact lines by the result for the fluctuation force per length / 2 d(^) between 
two parallel lines with a separation distance h and integrating over the two opposite contact 
lines to obtain the total effective force between the colloids. The Casimir force density between 
two parallel surfaces was calculated in Ref. 23( in a general approach for arbitrary dimensions. 
Applied to two dimensions we obtain the force line density f2d{h) = — k^T n / (24h 2 ) . Integrating 
this density over the opposing contact lines yields [24] 

■2 



24 



dy 



1 



-k-aT 



Ti- 



ll 



yif-L -l 1 

2 Vft T R 2 . 



48^1 V i + i 



0(h 



-1/21 



(22) 



Here, R\ and i? 2 are the curvature radii of the two ellipses at the end points of the distance 
vector of the closest approach. It is seen that the fluctuation force diverges as h~ 3 ^ 2 upon 
contact of the ellipsoids (h — > 0). 



A. Intermediate distances: Numerical calculation 



For intermediate distances d — d c \ ~ a the fluctuation induced force has to be calculated 
numerically. This can be done in principle by including a number of multipoles in the numeric 

11 



evaluation of the determinant in Eq. (Tl9|) . see Ref. [12|. In order to avoid the algebraic evaluation 
of the multipole coefficients of G, it is possible to apply a method which was introduced in 
Ref. [ijj]. The starting point is Eq. (fit)]) for the partition function Z m . Introducing a mesh 
with N points 77^ , j — 1 . . . N, on the reference contact line <9Si ire f converts the double integral 
in the exponent to a double sum. Thus the functional integrals over the auxiliary fields ipi{r]i) 
are replaced by ordinary Gaussian integrals over the ^(x^?]^)), Vipi ~ rijlo ^0i( x i (%'))■ m 
the exponent, the ^(x^^jj)) are coupled by a matrix G with elements G\i = G(\x.i(r]ij) — 
Xi'(Vi'j')\)- performing the Gaussian integrals and disregarding divergent and <i-independent 
terms immediately leads to jF in = {k^T/2) lndet(G^ 1 G(<i)) for the fluctuation free energy. 
Here, G^ = lim^oo G(d). It contains the self energy contributions and is needed for the 
regularization of the free energy. Deriving with respect to d, the Casimir force can be written 
as 

= -^tr [GW^dMd)] • (23) 

The advantage of the direct calculation of the force is that Eq. (!2B1 does not contain any 
divergent parts which would require regularization, thus easing the numerical treatment con- 



siderably. The determinant is computed by using a standard LU decomposition (25|. We find 
good convergence of the numerical routine. The convergence can be sped up by distributing 
more points in the regions where the ellipses face each other. We note that computing the force 
by the multipole series seems to be more efficient this can partially be compensated by 
the point distribution on the ellipses. In Fig. [3b, (ellipse aspect ratio 2) and [3b (aspect ratio 
6) we compare the analytical results of Eqs. (120!) and the Derjaguin approximation (Eq. (1221 ) 
with the numerical results. As it is shown the analytical expressions show very good agreement 
with the numerical data points for both long- and short range behavior and almost cover the 
whole distance regime. At large distances d, the leading term of the free energy expansion in 
Eqs. ( 120]) mainly determines the behavior of the Casimir interaction because of its long-ranged 
nature, hence the orientation dependence of the subleading terms can be neglected. In order 
to demonstrate the anisotropy of the Casimir interaction, we show results for a fixed, inter- 
mediate distance d between ellipsoid centers and varying orientation 9 2 of the second ellipsoid, 
see Fig. [3b (aspect ratio 2, d/b = 4.1) and(3U (aspect ratio 6, d/b = 12.1). The orientation of 
the first ellipse was fixed to three values, Q\ — 0, Q\ — 7r/4 and Q\ — tt/2. As can be seen, the 
fluctuation-induced interaction is maximally attractive for 9\ = 6 2 = (tip-to-tip configura- 
tion). When 62 deviates from zero then the resulting force reduces. This behavior holds for 
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FIG. 3: (a, b) Comparison of the numerical results for the interface fluctuation Casimir force (symbols) 
with the analytical expressions in the asymptotic ranges of large colloid separations (i> a,b (full line) 
and small surface-to-surface distance h = d — d c i <C & (dashed line) (c,d) numerical Casimir interaction 
between two fixed ellipsoids trapped at the interface as a function of their orientation, for d = 4.1 and 
d = 12.1, respectively. 

both aspect ratios 2 and 6. 



IV. INCLUSION OF COLLOID FLUCTUATIONS 



In general, the inclusion of colloid height and tilt fluctuations into the partition function 



(Eq. (jSj)) can be realized by an approach used in Refs. 0, 



lo| . In this approach, the partition 



function is split into a product of a colloid fluctuation part and the interface fluctuation part. 
The latter contains only the contribution of the fluctuating interface, with Dirichlet boundary 
conditions on the colloid surface (see previous section). In the colloid fluctuation part, the 

13 



fluctuations of colloid heights and tilts are weighted by a Boltzmann factor which contains 

the energy of the mean field solution (Euler-Lagrange equation) to the capillary problem with 

the boundary conditions set by the fluctuating contact line. This decomposition is possible 

due to the fact that capillary wave Hamiltonian is Gaussian in the field u. In principle, it 

is possible to use this method also for the special case of ellipsoidal colloids considered here. 

However, finding the mean field solution in such a geometry for arbitrary contact lines is rather 

cumbersome. To bypass this difficulty, we employed a trick adapted from Ref. [?J in which 

effective forces between rods on fluctuating membranes and films have been investigated. We 

extend the fluctuating interface height field u(x, y) which enters the functional integral for Z to 

the interior of the ellipses S^i- Thus the measure of the functional integral for Z is extended 

by XHt(x)| xg £. rcf and the integration domain in the capillary wave Hamiltonian is enlarged to 

encompass the whole R 2 . On the colloid surfaces, the interface height field is given by the three 

phase contact line, u(dS itTe f) = /j. We extend u continuously to the interior of the circles S ijTe f. 

Such a continuation is not unique. However, the partition function remains unchanged (up to 

a constant factor), if the energy cost of such a continuation is zero (as it is physically required 

since the interface is pinned to the ellipsoid surface). This has to be insured by appropriate 

counterterms 0, [k]]. We choose the continuation: 

in \ r 1+ \ \ ( r-, cos(/im£j) sinh(m£j) . . A 

u{S ijie {) = /j,cxt(4i, Vi) = ) A p im — f. \ cos(m?7;) + Q im sm^m^) , (24) 

\ cosh(m^ ) smh(m£o) / 

where ^ and r)i are the elliptic coordinates with respect to ellipse S^ref- The specific choice 
above is convenient for the further calculations since V 2 /i )ex t = in Si^ e {\dS itTC f. Extending 
the integration domain of the capillary wave Hamiltonian in Eq. dSJ), Q = R 2 \ UiSi tTe t — > 
R 2 generates an additional energy contribution — 7ij >C orr which has to be subtracted from the 
extended capillary wave Hamiltonian 7i cw [f2 = R 2 ]. Therefore, the total Hamiltonian reads: 

2 

Titot = 'Hew + y^[7Yj,b + ^i,corr] • (25) 
i=l 

The correction Hamiltonian is calculated in App. [B], and we recall the boundary Hamiltonian: 
-n hCorr = ^^m(^ 2 m tanh(me ) + QLcoth(me )) • (26) 

m 

n i>h = ^(f» 2 1 tanh(e )+Q? 1 coth(eo)) ■ 

In Eq. (1261) we have already omitted the contributions from the gravitational term in 7i cw which 
are of order (a/ A) 2 1. 
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As in the previous section the partition function is written as a functional integral over all 
possible configurations of the interface position u and the boundary lines, expressed by fa, 

Z = Z, 1 Jvufljvf, J] 5[u(^) - / iiext ( Xi )] exp {~ ^f^ } , (27) 

where the product over the 5-functions enforces the pinning of the interface at the positions of 
the colloids. In contrast to Eq. (jSj), this product extends over all x e S^ref instead of dSi iie {, 
only. The 5-functions can again be expressed by auxiliary fields ipi, now defined on the two- 
dimensional elliptical domains Si >ie { as opposed to the auxiliary fields of Sec. II III which are 
defined on the one-dimensional ellipses dSi Te f m . 



Z 



Vu J fl V ^j Vfiexp l~ 1 ^F^ +i J s ^^( x )Kx)-Ac Xt (x)]| .(28) 

Similarly to the evaluation of the fluctuation part, Sec. 1111} we introduce multipole moments 
ty im of the auxiliary fields by inserting unity into Z, Eq. ([28]) : 

1 = / II II ^i^L 5 (%n - I d * x (cosh(mO/ cosh(m£ )) cos(mr/) ^(x) 

x ^ l^im ~ / d 2 x (sinh(m£)/sinh(m£ )) sin(m?7) ipi(x.) J . (29) 

\ J <Si,ref / 

In contrast to the evaluation of the fluctuation term in Sec. 1111] there will be constraints on 
the lowest multipoles which contribute to Z. To see this we note that the Hamiltonian 7i tot 
does not depend on the boundary monopole moments P^ and the dipole moments Pn (through 
a cancellation between H.^ and Ti,i jCorv ), and the only dependence of Z on these moments is 
through the constraint function fi tCxt - Recalling the definition of the integration measure T>fi 
for the two boundary conditions (B) and (C) and performing the integration over P i0 (B) and 
Pio and Pn (C), we immediately find 



Z 



r 2 

J IIIl^^-W- case(B) 



i=l m 
2 



/ I! II d ^m ■ ■ ■ mo) mi) W • • • cases (C) 



(30) 



i=l m 

Having noticed these constraints on the auxiliary fields, we proceed by integrating over the field 
u in Eq. (p8l): 

2 

d 2 Xj ^j(Xj) CrflXj - X,-|) 1pj(*j) 



Z = [ f[v^ ( VU exp [ I 

J i= i J I 7 ■■ 1 Js i , tei Js 
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~7f (Hi.b + ^i.corr) ~ 1 / d' 

Bl i=l J s i>Te( 



a?'0i(x)/ i)e3rt (x) | , (31) 

where - as in Eq. (TTOl) - G is the Greens function of the capillary wave Hamiltonian. A 
somewhat longer calculation shows that Z can be split into into an interaction part (coupling 
the auxiliary multipole moments , P% m an d Qim for different colloid labels i), a self-energy 
part (depending on , Pj m and Qi m for each value of i separately) and a remainder (the sum 
of boundary and correction Hamiltonian): 

z = J nn^-^L/ ^ex P |-^(w int [*s,*£]+^seif[*s ] ])} x 



6XP I ( ^' b + ^' corr) " 1 E(*™ P - + I ( 32 ) 



III 



The interaction part 

Him =2 d 2 x x / ci 2 x 2 ^i(xi)G int (|xi - x 2 |)^2(x 2 ) (33) 

J Sl,ref J >S2,ref 

turns out to be a bilinear form in the auxiliary multipole moments; this is shown using the 
already used multipole expansion of the Greens function C(|xi — x 2 |) — — m(7 e |xi — x 2 |/2A c ) 
(valid for d ^> a) which is presented in App. Oin more detail. This bilinear form reads 

■ ■,. ( ) ( ^ ' ) ( ) (34) 

where the submatrices Gf^ t , Gfn t and Gg t have already been encountered in the calculation of 
the fluctuation part and are given by Eqs. (fT5l)- (fT8l) . The self-energy part (different from the 
corresponding one in the calculation of the fluctuation part) is evaluated in App. (Dj with the 
result 

ln( 7e (a + &)/8A c ) iTfC 2 
rti, se u — W i0 

1^1/ ^ c 2 l£ s 2 \ 
_|__ X / 31 im _|_ im j 

it £-^ o m \ 1 + tanh(m£ ) 1 + coth(m£ ) / 
Combining Eqs. fl32l) . (133]) . and (135]) . the partition function can be written as 
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(35) 




where the vectors M/j = (\&£ , Pm, VP?!, Pa, *& a , Qn, . . .) - in contrast to in Sec. [TTT] — contain 
all involved auxiliary and boundary multipole moments. The elements of the matrix H describe 
the coupling of these multipole moments, where the self-energy block couples multipoles defined 
on the same ellipses Si jTe { . Thus the diagonal part of the self energy matrix H se if can be read off 
Eq. (1261) and Eq. (|35|) while the off-diagonal part is determined by the term — i J2 m (^i m Pim + 
^imQim) in Eq. fl32l) . The elements of the interaction matrix H int are determined by the 
interaction energy 7Y int in Eqs. (1331) and (1341) and couple the auxiliary multipole moments of 
different colloids. All matrix elements representing couplings of other multipoles are zero. 
Similar as in Eq. ( fl2l) . the exponent in Eq. (J36l) is a bilinear form, however, here combined 
for all types, boundary multipole moments P im , Q im and auxiliary multipoles \l/f m , \l/| m . The 
computation of the partition function amounts to the calculation of det H. Again this is found 
as a series expansion in a/d, and we may define a similar expansion for the free energy ^in+coii = 
-(k B T) In det H/2 as before in Eq. (|20|): 

•Fin+coiiW = k B TY,fZ +c ° XX Q) • (37) 

The leading coefficients in case (B) (inclusion of fluctuations in the ellipsoids' vertical positions) 
are given by: 

/l n+co11 = ~ [ (a 2 - b 2 ) 2 005(29, + 26 2 ) + (a + bf ] . 

In case (C) (inclusion of fluctuations in the ellipsoids' vertical positions and tilt angles with 
respect to the interface) the leading coefficients are: 



fZ +C ° XX = (n = . . . 3) , (39) 
9 
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/8 n+C ° U = "4 [ (« 2 " b2 T cos ( 4 ^i + ^) + (a + b) 8 ] . 



In contrast to the calculation before, the different leading power laws for the different cases (B) 
and (C) can be understood easily. We note that the interaction between the auxiliary multipoles 
and m %nt: Eq. fl33l . scales like (a'/4d) m+n . After calculating the determinant, the 
leading order of the total fluctuation induced force between the two colloids is determined by 
the first non-vanishing auxiliary multipole moment and (as follows from det H) gives rise 
to a term in the free energy oc l/d im ' (for m' > 0) or cx lnlnci (for m' — 0). As explained in the 
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beginning of this subsection, the different boundary conditions lead to certain constraints on the 
auxiliary multipoles: According to Eq. fl30|) . the leading term in F(d) arises from a monopole- 
monopole interaction of the auxiliary field in case (A), from a dipole-dipole interaction in case 
(B), and from a quadrupole-quadrupole interaction in case (C). The constraints of vanishing 
auxiliary monopole and dipole moments (as in (C)) result from the independence of Htot of the 
boundary monopole and dipole moments and this is only captured correctly by the inclusion of 
the correction Hamiltonian 7i cor r- The Casimir attraction is maximal if the major axes of both 
ellipses are oriented parallel, regardless of the orientation of the distance vector joining their 
centers. This is a peculiarity in two dimensions, as can be seen also by the general multipole 
expansion of the interaction between two arbitrary charge distributions in two-dimensional 
electrostatics. 



A. Limiting cases 

In the limit a = b (colloids with circular contact line such as disks and spheres) our results 

nn 

for the cases (A)-(C) reduce to the results reported in Refs. |9|, [10[. In the limit b — > 
(colloidal rods or needles with vanishing thickness) we can compare our result for case (C) 
(fluctuating colloid heights and tilts) to Ref. Q. There it has been found that the effective 
free energy asymptotically varies oc d~ 4 with a coefficient given by Eq. fl38l) . i.e. by the result 
of case (B) (colloid height fluctuations only). The derivation in Ref. 3] suggests that the 
perturbative treatment employed in our approach should be amended by corrections in the 
integration measure over the tilts. In our cases, this measure is simply given by dPndQn, 
the product of the measures for the cosine and sine dipole moments of the contact line. In 
order to check the validity of this approximation for the measure, we recalculated the partition 
function of Ref. [3] by perfroming the integrations over the auxiliary dipole moments and their 
conjugate variables, which results in a final integral over the tilts (with the general measure) 
weighted by an exponential function in the tilts. If 7a 2 / '(k^T) ^> 1, the denominator of the 
measure in this integral can be expanded in terms of the tilt angles since the exponential 
function decays much faster than the measure and thus determines the convergence of the 
integral. In this way, one sees that the higher-order terms in the dipole tilt measure do not 
provide another leading behaviour in 1/d in the partition function compared to the leading 
quadrupole-quadrupole interaction which arises in our perturbative picture. A breakdown of 
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our perturbative treatment can be expected if the length of the rod a approaches the molecular 
length scale. This coincides with a simultaneous breakdown of the simple capillary wave picture 
underlying our analysis. 

V. SUMMARY 

The restrictions that two colloids trapped at a fluid interface impose on the thermally excited 
interfacial fluctuations (capillary waves) by their sheer presence lead to a thermal Casimir 
interaction. We have obtained an explicit account for the effect of colloidal anisotropy on 
the form of the Casimir interaction by studying ellipsoidal (spheroidal) colloids with arbitrary 
aspect ratio. For the case of fixed colloids and fixed contact lines, the problem is equivalent to 
the "standard" Casimir problem for a scalar, Gaussian field in two dimensions with Dirichlet 
boundary conditions on the colloid surface. In an expansion in 1/d, the inverse center-to-center 
distance between the colloids, the leading term in the Casimir interaction energy is found 
to be attractive, isotropic in the interface plane and slowly varying oc lnlnci (see Eq. (EE))). 
Anisotropies appear in higher orders in 1/d and become important when the closest surface- 
to-surface distance between the colloids becomes small (see Fig. [3]). 

If fluctuations in the colloids' vertical position are permitted, the asymptotics of the Casimir 
interaction energy changes to a behaviour oc d~ 4 (see Eq. ( |38l) ). In this case, anisotropies 
are present in the leading term but the interaction remains attractive for all orientations. If 
furthermore fluctuations of the colloids' orientation with respect to the interface normal are 
allowed, the asymptotics changes to a behaviour oc d~ 8 (see Eq. ( 1391) ). Interestingly, this 
change of leading order in the asymptotics of the Casimir energy depending on the type of 
permitted colloid fluctuations holds for arbitrary aspect ratios. This leads to the speculation 
that this might be a general feature holding for arbitrary colloid shape. 

In our approach, the Casimir interaction can be understood as the interaction between fluc- 
tuating multipole moments of an auxiliary charge density-like field defined on the area enclosed 
by the contact lines. These fluctuations are coupled to fluctuations of multipole moments of 
the contact line position which are a due to the possibly fluctuating colloid height and tilts. 
Therefore, the system can be viewed as an example for the Casimir effect with fluctuating 
boundary conditions. Such fluctuating boundary conditions appear to be difficult to be real- 
izable in three-dimensional systems such as the standard system of charged metallic objects 
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subjected to vacuum fluctuations of the electromagnetic field. 

Experimentally, the detection of the Casimir interaction at a fluid interface appears to be 
possible if competing interactions, especially van-der-Waals and static capillary interactions, 
are sufficiently weakened. Van-der-Waals interactions are also strongly attractive at small 
distances, but can be modified by an appropriate core-shell structure of the colloids or by 
using flat, disk-like particles (we refer to a longer discussion of this issue in Ref. 9[). Capillary 
interactions are very strong for ellipsoidal colloids of micrometer size and with contact angle 
different from 7r/2 since the equilibrium contact line in this case is already undulated and gives 



rise to static deformations of the surrounding interface [24j]. These capillary interactions can 
be minimized by either using truly nanoscopic ellipsoids or synthesizing Janus particles with a 
contact line which is flat on a nm level. Despite the great advances in particle synthesis over 
the last years, this appears to be still a big challenge. 

Acknowledgment: The authors thank the German Science Foundation for financial support 
through the Collaborative Research Centre (SFB-TR6) "Colloids in External Fields", project 
N01. 

APPENDIX A: CONFOCAL ELLIPTIC COORDINATE SYSTEM 

Confocal elliptic coordinates (£, rj) are planar orthogonal coordinates formed by confocal ellipses or 
hyperbolae. The foci are located on the x-axis of the Cartesian coordinates, separated by a'. The 
relation to Cartesian coordinates is defined by 

x = — cosh(^) cos (77) , 

V = % sinh(£) sin(r?) , (Al) 

and the scale factors are found as 



a' 2 

h^ = h v = d — (cosh(2£) - cos(2? ? )) . (A2) 

£ and r\ are called elliptic radius and elliptic angle, respectively. In this coordinate system, £ = £0 

represents the equation of an ellipse with axes a, b (a > b). The elliptic radius and the distance a' 

between the foci are given in terms of the ellipse principal axes by 

, 1. (a + b 
£0 = -x ln 



2 \a-b / 

a' = (a 2 -6 2 )t. (A3) 
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Therefore, for a circle we have £o — > 00 an d for a line, £0 ~ 0. 



APPENDIX B: THE BOUNDARY AND CORRECTION HAMILTONIANS H i)h AND 

^i,corr 

i) The boundary Hamiltonain for the case of a pinned contact line (Janus ellipsoids) is governed 
by the difference in projected meniscus area, Aj4 pro j (Eq. ([2]). If the ellipsoid is tilted in the 
xz-plane by a small angle oti, this area is given by 



A^w— a\a sinh(2£ ). (Bl) 



The contact line position u\gs- T { is easily determined by using a coordinate system rotation 
from the xz-plane to the x z -plane, where the x - and the z -axis coincide with the major and 
minor axis of the tilted ellipsoid. The contact line is located at z = in the new coordinate 
system. Since z ~ z — ajX, we find that u|g5 iref = z\ z '=a ~ a i x an d therefore u\gs iTcf ~ 
aj(a'/2) cosh(^o) cos(?7i). Thus, in the multipole expansion of the tilted contact line (Eq. ©), 
only a dipole term appears with the dipole moment given by: 

P a = ai yCosh(Co) • (B2) 



Inserting eq. (|B2[) into (jBTJ) we obtain 

A^ oj = |p 2 tanh(£o). (B3) 

Applying the same arguments for tilts in the yz-plane, we can express the boundary Hamiltonian 
in the small-tilt approximation by 

%,b = 7AA pro j 

= ^{Pl tanh(£ ) + Ql coth(Co)) . (B4) 

ii) The correction Hamiltonian, which is introduced in Eq. (|25p . is determined by minus the surface 
energy of the meniscus piece u\s iiel = f (Eq. ([24]) ) extended into the ellipses enclosed by the 
reference contact lines: 

,,2- 



7 / 2 



(Vuy + ^ 



Ac =°° I _/ d 2 ^(V/) 2 . (B5) 
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Thus, 

f<0 /-27T 



Wi,«*r = 1 I" / 7r ^r ?i ((%/) 2 + (^ i /) 2 ) 
z Jo Jo 



^^m(^ 2 m tanh(m^ )+Q 2 m coth(mCo)) . (B6) 
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APPENDIX C: EXPANSION OF GREEN'S FUNCTION IN ELLIPTIC COORDI- 
NATES 

In this appendix we derive the multipole expansion of the Green's function G(|x|) ~ 
— (l/27r) ln(7e|x|/2A c ) between two charged elliptic regions (charges are generated by auxilliary fields 
ipi). This Green's function gives the correlation between two points residing either on the same ellipse 
or different ellipses. 

1) In the case that xi and X2 are located on the same ellipse, the Green's function expansion has been 



given in Ref. 



29| : 



27rG(|xi— X2I) = — In +2 y [cosh(n^o) cos(nr/i) cos(n?72) + sinh(n£o) s'm(iir]i) sinfra^)] 

(Cl) 

2) The case that the two points are located on different ellipses, i.e. xi = ri and X2 = d + r 2 , which 
furthermore possess an arbitrary orientation in the plane (expressed by the angles 8\ and 82, see Fig.Q]) 
is more difficult. We start with a general Taylor expansion of the Green's function: 



2vr V 2A C J 27T \2\J 27T j,\ j 2 \ 



• (C2) 

r=d 



On the other hand, we can perform this expansion using complex variables z = x + ly. The expansion 
of the logarithm is given by: 

00 ^ 

hx(z-z) = Y -(-z'd z ) j In z . (C3) 
j=o J - 

The real part of eq. (|C3h is the expansion of the real logarithm, of course: 

ln|r-r'|=Ref;-i^. (C4) 

3=0 J 1 1 

By comparing Eq. (|C4p and the Taylor expansion of In |r — r | in real space as in Eq. (|C2p . we find 

(zi_ZH lllr = _I Re (^)i. (C5 ) 
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Introducing complex derivative operators (_ = d z , = d z * and using C+C- lnr = C-C+ lnr = 0, we 
have 



Cilnr = Ciln|z| = ^(-l)^ 1 (i-l)! < 



^7kl 2j 



(C6) 



Identifying eq. (|C5|) and (|C6p we obtain 



^LZ)ilnr = -i ( (z *) J C1 + P'Cl ) lnr . (C7) 



By inserting eq. (|C7p into eq. (|02|) . we find the Green's function expansion in terms of complex 
variables Zi 

1 , ficd\ l (- 1 )- 72 ( h+h\ 1 „ r ,- 91 

G x 2 -xi = In ^- + — > ^ — '— -—. Reh^d 2 ] . C8 

31 + 32 > 1 

This general expansion can be used in any coordinate system. In the special case of elliptic coordinates, 
z = (a /2)e~ %e cosh(£ + i?7), where 6 is the in-plane rotation angle of the ellipse major axis with respect 
to a fixed x-axis. (For the configuration of arbitrarily oriented ellipses, the x-axis is given by the line 
joining their centers, see Fig. [I]). In order to express Eq. (|C8j) in elliptic coordinates, one needs z 3 
which can be found by applying the binomial expansion: 

■ _ / a e~ %9 \ A / j\ (2k _ m+iv) 



k=0 



[ J e (2fc-i)(€+*j) . (C9) 



Using the above expansion, Eq. (|CJ8|) becomes: 

c( , XI _ X2 | ) = __L ln (g 



+^ E EE 



(-1) J2 (h +h\ ( ji\ ( h\ [ a 



27F i in i n + J' 2 V h J \klj \k 2 J \4d I 

31,32— V fel=0fc2=0 \ / 

31+32> 1 

x exp(ei(2/ci -h)+ 6(2*2 - J'2)) 

x cosOi^! + j 2 2 + % (ii - 2fci) + r? 2 (j 2 - 2fc 2 )) . (CIO) 

The aim is to rewrite this fourfold sum over j\ , j'2 , k\ , fc 2 as an expansion into multipole coefficients 
cos (777,771) cos (71772) and sin(m?7i) sin(n772) with m, n > 0. To that end, we define m = \ji — 2k±\ and 
n = \]2 — 2Jt2 1. The possibility that jj — 2ki may be positive as well as negative makes it necessary 
to consider the following cases: (i) "m = 0, n = 0", (ii)"m = 0, 0" or "m 7^ 0, n = 0" and (iii) 
u m 0, n / 0". 
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(i) "m = 0, n = 0": Two auxiliary variables l\, Z 2 are introduced through j% = 2k,ki = k (i = 1,2). 
The such constrained sum in Eq. ([ClOp reduces to 

£^r* + X)£)(C"^ +sw - 

i 1 +i 2 >0 

Relabelling I = l\ + I2 and I' = h, the above sum is rewritten as 

1 1 
E E 97 



©(T)(1-") (a) <-> 



(ii) "m = 0, n > 0": Here, Zi is introduced as above through j\ = 2l\ and k\ = l\. We distinguish 
the two cases ji ~ 2^2 > and j 2 — 2ki < via the choice of I2 through ji = n + 2Z 2 and fc 2 = Z2 
vs. &2 = n + I2. Adding up these two cases in the constrained sum (|C10p . and performing a 



relabelling analogous to the one leading to expression (jClip (I = h + h, I' = h) we obtain: 
r(n + 2/) / a' 

'/' 2 (Z-Z')!(n + Z-i')! ^ 



2 / ' \ 

2 E E(- 1 ) n 7^77ZW7 : T7377T7 h {cos(2l'e 1 + (n + 2(1 - l'))9 2 ) cosh« 2 ) cos(n % ) 



+ sin(2Z'6'i + (n + 2(1 - l'))6 2 ) sinh« 2 ) sin(nr? 2 )} . 

(C12) 



Similarly we obtain for "m > 0, n = ": 

jr / , \ m+2l 

2 E E - fOl^mV ZO! I {C ° S((m + 2/>1 + " °' 2) COSh(m6) COS(m?? 



+ sin((m + 2l')6i + 2(1 - l')0 2 ) smh(m£i) sin(mf?i)} . 

(C13) 

(iii) "m > 0, n > 0": Similarly to the previous cases, ji and j 2 are introduced through ji = m + 2l\ 
and j 2 = n + 2Z 2 . The four cases of possible sign combinations of j\ — (i = 1, 2) are taken 
into account by the relation sets tl k\ = l\ , & 2 = ^2", " &i = h , k 2 = n + Z 2 ", "fci = m + Zi , fc 2 = 
fa""ki = m + l\, k 2 = n + Z 2 ". Adding up these four cases in the constrained sum (jClOp and 
taking advantage of the addition and subtraction relations between hyperbolic functions and 
then relabelling as before (I = l\ + Z 2 , V = l\) we find: 

4 V Vf-lV 1 r(m + n + 2Z) / o_\ 

(m + ZOlZ'Kn + Z-Z')!^-^)! \ 4 V 
(cos(G) cosh(m^i) cos(mr7i) cosh(ra£ 2 ) cos(nry 2 ) 

+ sin(G) sinh(m£i) sin(mryi) sinh(n£ 2 ) sin(nr/ 2 )) , (C14) 
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where G = (m + 2l')6 1 + (n + 2(1 - 0)^2- 



The cases considered above may be combined into a single expression such that the 2D Green's function 
in elliptic coordinates reads: 

/ / \ m+n+2l 
v 7 m,n,l=0 \ / 

x {A mn i [cos(mr/i) cosh(m£i) cos(nr/2) cosh(n^2) — sin(mr/i) sinh(m£i) sin(n?72) sinh(n^2)] 
+A s mnl [sin(m77i) sinh(m£i) cos(nry2) cosh(n^2) + cos(mryi) cosh(m£i) sin(nry2) sinh(n^2)]} • 
In Eq. (IC15j) . the coefficients are given by 

^ = AB ll (0u0 2 ) = 2^pcr , (C16) 

where Hj is the discrete step function; H j = 1 — 5jo, and 



rimn _ i i\n T(m + n + 2/) / P17 i 

C "' " (m + l')W\(n + l-l')\(l-l')\ ' (L17) 



with C °£ = 0. 



APPENDIX D: SELF ENERGY IN THE CASE OF FLUCTUATING COLLOIDS 

For fluctuating colloids, the self-energy part of the multipole-multipole interaction of the auxiliary 
fields was introduced in Eq. f|32|) . Here we calculate it explicitly with a method similarly to the 
one employed in Ref. [3]. As the starting point, we obtain from Eqs. ([HTj) and (|52"|) the following 
expession for Z i>self = exp (-%^Wj, S cif) : 

x 5 (v s im - jf d 2 x (sinh(?0/sinh(m^o)) sin(mr ? )^(x)^ exp ^i^(*? m P im + *f m Q im )j 
xexpf-^/ d 2 x / dV^(x)G(|x-x'|)Vi(x)-i / d 2 x ^(x) / ij(Dtt (x) ] . 

\ 2 7 ./s iiref / 

(Dl) 

The (^-functions in 2j jSe if may be eliminated by introducing conjugate multipole moments \l/? m and 
Vl/f m to the multipoles ^ of the auxiliary fields: 

5 (v® - f d 2 x ( c °M^i)cos(m Vi )/cosh{mCo)\ \ = 

\ m Jsi, iet V smh ( ra ^) sin ( ?7i ^)/sinh(?n^o)y * ) 
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/ 



C) 



2 /cosh(m£;) cos(mr]i)/ cosh(m£ )\ 
d x 1 ■ w t • w • w m m x i 



(D2) 



Inserting Eq. (jD2[) into 2j se if we obtain: 

k B T 



-Zj.self = / JJ^im / V^exp 

cosh(m^j) 



d 2 x / d 2 xVi(xi)G(|x - x |)^i(xj 



il d x ipi(xi) 



cosh(m^o) 



„ ~ r . , . sinh(m^) ~ . 

-Pjm + Vw cos{mr]i) + -— - — ~-r(Q im + Vim) sm{mr]i 



sinh(m^o 



(D3) 



+i Yl + V C imWim + (Qim + ^im)% 

where d^j m = d^> s im d^ c im . The functional integral J T>^i in eq. (|D3h can be replaced by a functional 
integral over a constrained height field h(pc): 

-Zj.aelf = / II eXP ( 1 [( Pim + + (Qim + *!J*L] ) 

m \ m=0 / 

x Jvh \[ 8(h(xi)-fi)exp 



7 

2/c B T 



d 2 x 



(V/i) 2 + 



(D4) 



cosh(m&) ( , cosfmr/i) + sm M m ^) (q + j sin ( mJ7 

cosh(?n£ ) m smh(m£o) 



where /, = J2 r , 

In the region Si jTe {, i.e. the ellipse enclosed by the reference contact line, the height field h is pinned 
to /j. Therefore the contribution of the functional integral J T>h in this region is simply given by the 
surface energy of fi which was determined in Eq. (jB6[) and itself becomes: 

Z iM i Ac S°° J JJ d% im exp | ^ ™ [(^im + ^ m ) 2 tanh(m£ ) + (Qim + ^tmf eoth(m£ ) 



+ij^ [(P; m + VtmWim + (Qim + 
m 

x Jvh 11 5(/ l (x i )-/i)exp ' 



d 2 x 



h 2 



(D5) 



where the remaining 5-functions describe the pinning of h(x) to the boundaries dSi ze t of the integration 
domain. The auxiliary field can be separated into two parts, h = ho + hi, where (—V + A~ 2 )/io = 
with the boundary conditions ho(xi)\gs i fef = fi and /ii(xj)|a5 i rcf = 0. Applying Gauss' theorem to 
the integral in the exponent of Eq. (|D5|) leads to 

/ ]J eXP ( ~2k^T 2 + tanh ( m ^°) + + ^irn) 2 COth(me ) 



m 

x j Vh x Yl 6(hi(xi))exp 



exp 



7 



2/c B T 



dx. /io(x)V/iq(x) 



2/c B T y R 2\ S . ref 



(f 2 x 



(v/m) 2 + | 



(D6) 
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The functional integral over h\ yield to a constant value independent of any multipole moment, which 
can be neglected. 

The general solution to the Helmholtz differential equation for ho in R 2 \ Si jre f is needed for computing 
the line integral in eq. (|D6|) . Separation of variables in the Helmholz equation in elliptic coordinates 
would lead to the Mathieu (angular part) and modified Mathieu (radial part) differential equations. 
The solution to these equations in the asymptotic case A c — > oo become standard triangular functions 
(angular part) and modified Bessel functions of the second kind (radial part), respectively. Therefore, 
the solution is given by 

ho(xi) = Y] ^^> e J?n\\ ( Am cos ( m Vi) + B m sm{mr]i)) . (D7) 
^ K m (a e*o/2A c ) 

The coefficients in Eq. (|D7j) are readily determined by comparing to the boundary conditions: A m = 
Pim + ^fmi B m = Qi m + Vl/| m and Bq = 0. Then the line integral evaluates to 2Trg(m)(A^ n + B^) with 
g(m) = m/2 (m > 0) and g(0) = — 1/ ln(7 c a'e^°/8A c ). Using this, ^ se if reads: 



i.self 



/ ]J d ^ im exp ( ~uTT 2 m [ (Pim + ^ m)2 tanh ( m &) + ( Qim + ^imf coth(m£ ) 

+i^ [(p im + njn* + (Qim + *i m )*i m ] ) (D8) 



in 



x exp \~1Tt ^ 5(m) [ (Pim + ^ m)2 + ( Qim + 

V B m=0 / 

The last integration over the conjugate multipole moments can be performed after shifting variables, 
~ *■ + ^im an d *fm — * + *im" After this integration, the final result for H^ Be i{ = 
(-2j)/(k B T) lnZ ijSclf is given by: 

rti,self = 

_|__ _ j ™. | Y im \ (D9) 

7r m \1 + tanh(m^o) 1 + coth(m£o) / 
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